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Abstract 



In this chapter we explain the principles of velocity estimation and review both analytically and using 
simulations the properties of existing estimators. We then discuss the importance of accurate velocity 
estimation in the context of handover algorithm design. We show how an error in velocity estimation can 
significantly increase the probability of dropped calls, leading to a poorer quality of service of the system. 
Some velocity estimators require prior estimation of the Rician factor. This chapter also presents two Rice 
factor estimators that have been shown to have better performance than existing estimators. 
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5.1 Introduction 



This chapter presents the principles of mobile velocity estimation and illustrates its importance in the 
context of handoff algorithm design for mobile communications systems. Current mobile communications 
systems are under constant pressure to increase their capacity while maintaining a high quality of service. 
This is due to an ever-growing population of mobile phone users and an increasing demand for multimedia 
services. The main constraints system designers have to face are the distortions introduced by the mobile 
communications channel. These can be classified as Rayleigh or Rician fading, shadowing, and path loss 
[62, p. 16]. 1 These distortions limit the performance of mobile communications systems in terms of 
efficiency and quality of service. Many standard functions of the system are significantly enhanced by a 
priori knowledge or estimation of the mobile velocity. The next section provides examples where mobile 
velocity plays a key role in the overall performance of the communications system. 

5.1.1 Importance of Velocity Estimation 

In microcellular systems, the cell size is smaller than in macrocellular systems (Figure 5.1). As a result, 
microcellular systems require faster and more reliable handovers. 2 Also, since the base station (BS) is at 
lamppost level, if the mobile station (MS) rounds a corner, the signal power received by an MS can drop 
rapidly by 20 to 30 dB over distances as small as 10 m (Figure 5.2 and Figure 5.3) [6]. In that case, an 
emergency handover needs to be processed toward a target BS if the call is to be maintained. The drop 
in signal strength can be detected by applying short temporal window averaging in the received signal. 
However, the window size is velocity dependent and is optimal only when an accurate estimation of the 
mobile velocity is available [7]. 

In multitier systems (Figure 5.1), cells of different sizes coexist in a two-layer structure, i.e., microcells 
on the lower layer and macrocells (umbrella) in the upper layer. Within these systems, different types of 
handover have to be managed between the umbrellas and the microcells. In order to minimize the number 
of handovers, the MS velocity can be used in a cell layer assignment strategy, in which an MS is allocated to 
different hierarchical layers according to its velocity. The umbrella cells are used for fast-moving users and 



'The distortions of the mobile communications channel are reviewed in Section 5.2. 

2 Handover or handoff is the process of transferring the control of a mobile station from one base station or channel 
to another. It is an essential component of mobile communications systems. 




FIGURE 5.1 A multitier (microcell/macrocell overlay) system. 
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FIGURE 5.2 Corner effect: the MS turns round a corner, the LOS front the current BS (BS1) is lost, and an LOS is 
established between the MS and the target BS (BS2). 



microcells for slow-moving users. This results in a reasonable grade of service (GoS) for both microcellular 
and multitier systems [33]. 

Reliable estimates of the MS velocity are also useful for effective dynamic channel assignment and the 
optimization of adaptive multiple-access wireless receivers [26, 69] . Since the performance of many receiver 
techniques depends on the fading rate of the received signal, an adaptive communications receiver can im- 
prove the performance and reduce the complexity of current systems by using Doppler information to con- 
trol the receiver parameters. These parameters include the pilot filter bandwidth, the automatic gain control 




FIGURE 5.3 Received signal strength (RSS) at the MS from two neighboring BSs as a function of distance. 
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loop bandwidth, the phase tracker bandwidth, and the size of the interleaver [ 19, 52] . Based on this idea, Lee 
and Cho [42] proposed a power control scheme for codedivision multiple-access (CDMA) systems, which 
selects the power control step size based on the MS velocity. The MS velocity also affects the performance 
of a communications system operating in a pseudonoise (PN) tracking in CDMA systems [20, 22], 

Mobile velocity estimation can occur at both BS and MS sites. For example, in the current CDMA 
2000 system, the BS is responsible for power control [4, 41], and velocity is therefore estimated by the 
BS. However, a decision to hand off is made by the MS, and the MS is then responsible for estimating 
the velocity. Velocity estimation is performed by exploiting the statistics of a received pilot signal. For 
example, in CDMA 2000 the pilot signal consists of a constant-amplitude, zero-phase signal, spread by an 
all-one code sequence. It is then added to the information signal spread by an orthogonal sequence. The 
combination of the two signals is modulated by a carrier frequency and transmitted [21]. This process is 
equivalent to transmitting a pure carrier of constant amplitude alongside the information signal. At the 
reception, the combination of the two signals is demodulated, then correlated by each spreading sequence. 
Assuming no loss of orthogonality, the correlation of the received signal with the all-one sequence cancels 
the information signal and reveals the distortions introduced by the channel on the pilot signal. 

5.1.2 Existing Velocity Estimators 

Several mobile velocity estimators have been proposed in the literature, and some have been reviewed in 
[65]. Among the existing estimators are the zero-crossing rate of the in-phase or quadrature component 
[ 7] , level crossing rate [ 7] , and the autocovariance of the envelope of the received signal [5]. In [30, 31], the 
received signal samples are used to estimate the autocovariance function of the received faded envelope, 
from which the velocity information is extracted. Similarly, in [61], two methods for approximating 
the mobile velocity are proposed, based on the autocovariance function of the envelope and quadrature 
components of the received signal. The method proposed in [55] is based on the squared deviations of the 
logarithmically compressed signal envelope. In [ 1 ] , the rate of maxima of the envelope of the received signal 
is used as a velocity estimator. The method in [49] is based on the estimator proposed in [ 1] and applies 
the continuous wavelet transform to locate the extrema of the received signal envelope. In [66], velocity 
estimators are derived, based on the spectral moments of the received signal. The method described in 
[36] uses the switching rate of diversity branches in a selection diversity combiner to estimate the MS 
velocity. In [68], mobile velocity is estimated by applying an eigenmatrix pencil method to the received 
signal samples. In [71], two methods are proposed for estimating mobile velocities for Global System 
for Mobile Communications (GSM) radios. The methods are based on estimating the deviation of the 
received signal strength. The method proposed in [50] uses the local stationarity of the received signal and 
expands it on a basis of smooth local complex exponentials. Velocity estimates have also been obtained 
by estimating the maximum Doppler frequency using eigenspace methods [6]. The method is designed 
under specific assumptions of the angular distribution of the incident power. Other velocity estimators 
have been proposed that require prior estimation of the channel [40] and covariance function of the 
channel power [48]. In [28], multiple BS and multidimensional scaling are used to estimate the velocity. 
This method requires knowledge of the average signal strength for all locations. In [26], a maximum 
likelihood estimator is derived and requires prior estimation of the channel parameters. A few estimators 
have also been proposed, which do not provide an explicit formula for the MS velocity, but rather classify 
the MS velocity as being fast, medium, or slow [39, 73, 74] . 

It was shown in [ 10, 1 1 ] that since all the above-mentioned techniques are based on the statistics of either 
the envelope or quadrature components of the received signal, they are not robust to shadowing. Other 
estimators have been proposed in [8, 9] based on the instantaneous frequency (IF) of the received signal. 
Recent work has already exploited the concept of IF estimation and time-frequency signal processing 
in wireless communications, e.g., in blind source separation, channel coding and capacity, interference 
excision, multiuser detection, code design, multicarrier transmission, and synchronization [2, 3, 12, 13, 
17, 27, 35, 37, 38, 56, 57, 58, 59, 60, 72], and led to significant improvement in the systems [16, Chapter 13], 
As shown in [9], the main advantage of IF-based velocity estimators is that they are robust to shadowing. 
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5.1.3 Structure of the Chapter 

The chapter is structured as follows. Section 5.2 focuses on the received signal model in a mobile commu- 
nications environment. In Section 5.3, we explain the principles behind velocity estimation and provide 
examples selected from the above-mentioned list of velocity estimators. The examples are selected on a 
basis of simplicity and performance. Section 5.4 provides the framework for analyzing the performance 
of velocity estimators with application to one specific estimator. In Section 5.5, we present simulation 
results of a selection of velocity estimators. Because some velocity estimators require prior knowledge of 
the Rice factor, we derive in Section 5.6 two Rice factor estimators. Section 5.7 discusses the importance of 
accurate velocity estimation in the context of handover decision algorithms, and Section 5.8 concludes this 
chapter. 

5.2 Received Signal Model and Statistics 



5.2.1 Received Signal Model 

Propagation, in a mobile communications environment, is subject to three phenomena: path loss, shad- 
owing, and multipath fading [62, p. 16]. Path loss affects all radio communications and refers to the loss 
of received power as the MS travels away from the BS. Path loss is completely characterized by the distance 
between the BS and the MS, the operating wavelength, the antenna height, and the surrounding terrain 
[62, p. 16]. Shadowing, or shadow fading, is caused by terrain configurations between the BS and MS and 
is produced by variations of the average of the envelope of the received signal over a few wavelengths. It 
is a random process and follows a lognormal distribution, i.e., when expressed in decibel, its distribution 
is Gaussian [62, p. 87]. Multipath fading, also referred to as fast fading, is due to the constructive and 
destructive superposition of many reflected, scattered, and diffracted plane waves arriving at the MS with 
different time delays and phase shifts. Figure 5.4 represents a typical received pilot signal envelope. 




FIGURE 5.4 Typical received signal strength at an MS traveling at v = 50 km/h. 



Copyright © 2005 by CRC Press LLC 





We assume, in what follows, that the distortions introduced by the channel are described by the model 
proposed by Lee and Yeh [44] . In this model, the real signal y(t) received by an MS (or a BS) is the product 
of path loss, shadowing, and multipath fading with some additive noise: 

y(t) = m(t)s(t) + n(t) 

= p(t)m 0 (t)s{t) + n{t) (5.1) 



where s{t) represents the multipath fading, m 0 (t) is the lognormal shadowing given by [49], 



m 0 (t) = 10 iW/2 ° 

L (f) is a zero-mean Gaussian process, and p{t) represents the path loss given by [49] : 

p(t) = P 0 (vt + d 0 )~ a/2 



(5.2) 

(5.3) 



Po accounts for antenna parameters, transmitted power, and other relevant system parameters; vf repre- 
sents the distance between the BS and an MS traveling at velocity v at time t ; d 0 is the distance between the 
BS and MS at time t = 0; and a, the exponent of the distance dependence, reflects the amount of power 
loss as a function of distance. The additive noise, n ( t ), is assumed bandpass Gaussian with one-sided power 
spectral density (PSD) S^if) given by [7]: 



S N (f) 



f : l/-/cl<f 
0 : l/-/cl>f 



(5.4) 



where B 0 represents the system bandwidth and is chosen equal to , where v max is the maximum mobile 

velocity and X is the operating wavelength. represents the maximum expected Doppler frequency 
over the range of velocities. 



5.2.2 Multipath Component Model 

The multipath component, s ( f ) , is a result of the superposition of a number N of incoming waves, including 
a possible line-of-sight (LOS) component. Each incident wave w k {t) is modeled as a constant- amplitude, 
time-varying phase signal: 

w k (t) = a k e’ l2 *f' t+ * k{t)) (5.5) 



where a k represents the amplitude of wave k and $> k (t) is the phase shift caused by the moving vehicle. If 
w k (t) arrives with angle of incidence 6 k with respect to the direction of travel of the MS, the Doppler shift 
created is 

fk = 7 cos(9 k ) (5.6) 

A 

and the subsequent phase shift introduced by that wave is 

<Mt) = 27tf k t + 'l > k = 2jT^-cosidk)t + T'r- (5.7) 

A 

If Wk{t) is a LOS component, 9 k and are deterministic and depend on the position of the vehicle. 
Otherwise, 9 k and are random and are generally modeled as uniformly distributed between (0, 2 jt) 
[34]. Combining Equation 5.5 and Equation 5.6, the multipath fading, s(t), can be written as 



;(t) =9t < 5> tW 



. k= 1 



= 9t < e j27tf ‘‘ 



E 

_k = 1 



a k e 



jljt l cosWk^+Vk 



(5.8) 

(5.9) 
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where 9i{}refers to the real-part operator, Wi(f) is a possible LOS component, and{w 2 (f), w 3 (f) . . . ,%(()} 
are the result of N — 1 independent scatterers and, as a consequence, are assumed independent and 
identically distributed (i.i.d). Signal s (f) can be further written as 

s(f) = S;(f) cos2jr/ c f — s q (t) sin2jr/ c f (5.10) 

where s,- (f ) and s q (f ) are the respective in-phase and quadrature phase components of s (f ), defined as 

N 

S{(t) ~ ak cos (lx- cos(6k)t + 4^) = Xi(t) + iiij (5.11) 

k= 1 
N 

s q {t) = '^2a k sin (in^cosie^t + = x q (t) + m q (5.12) 

k = 1 

where E {x,(f)} = E[x q (t)} = 0, hi,- = £{s;(f)}, m q = E{s q (t)}, and £{} represents the expectation 
operator. 

In the presence of a LOS component, the means m, and m q are nonzero [34] and 

s(f) = (xj(t) + nii) cos 2nf c t — (x q (t) + m q ) sin2jr/ c f (5.13) 

= x(t ) + m; cos2jr/ c f — m q sin 2 nf c t (5.14) 

where x,(f) and x q (t) are the respective in-phase and quadrature phase components of signal x(t), i.e., 

x(t) = Xi{t) cos2nf c t — x q (t) sin2jrf c t (5.15) 

5.2.3 The Scattering Distribution 

An important parameter in mobile communications system design is the distribution of incoming waves 
around a particular MS. It is referred to as scattering distribution. For example, in a typical macrocellular 
environment, the MS is usually uniformly surrounded by local scatterers, so that the plane waves arrive 
from many directions without an LOS component. As a consequence, the scattering distribution is usually 
considered isotropic. However, in a microcellular environment, the antennas of the BSS are only moderately 
elevated above the local scatterers. As a result, an LOS component may or may not exist and the scattering 
distribution is usually nonisotropic. The von Mises density provides a model for the distribution of the 
angle of arrival of the incoming waves [ 1 , 45 ] . It is a function of one parameter, x > which determines the 
directivity of the incoming waves. It is given by 

^ — T7 — i e/C ° S "> *>0, -tt<0<7T (5.16) 

2ttI 0 (x) 

where /„(.) is the modified Bessel function of order n and# is the angle of incidence of the incoming waves. 
Figure 5.5 shows the polar plots of p(9) against the angle of arrival of the plane waves for three different 
values of x ■ It can be seen that for x = 0 the scattering distribution is isotropic and that it becomes more 
directive as x increases. 

5.2.4 Statistics of the Multipath Fading 

When the number of incoming waves is sufficiently large (generally greater than six, [34, p. 69]), the 
in-phase and quadrature components of x(t), defined in Equations 5.11 and 5.12, respectively, tend to 
be independent zero-mean Gaussian processes, with variance a 2 . As a consequence, the envelope of s (f), 
defined as 
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|s(f)| = y/ (x ; (f) + mi ) 2 + {x q (t) + m q ) 2 




(a) x = 0 (b) x = 1 (c) % = 1 0 



FIGURE 5.5 Polar plots of p(9) in terms of 9 for x = 0, x = 1> and x = 10. For x = 0 the scattering distribution 
is isotropic and becomes more directive as / increases. 



follows a Ricean distribution with Rice factor K = 2L, where;; 2 = mf + m 2 q [62, p. 518]. The Rice factor 
K represents the ratio of the power in the LOS component and scatter components of the signal s(t). 

The PSD of x(t ) is given in [62, p. 42] as 

Sx(f> = ( VW<f^E !lGmpm + G(_9) r (_<>)1 ^ ' f-a<f- (5 . 17) 

[ o : \f - fc\ > fm 

where /„, = ^ is the maximum Doppler frequency shift, p(9) is the scattering distribution, and G(6) is 
the gain of the MS antenna. We will assume in the remainder of the chapter that a vertical mono-pole 
antenna with G(9) = | is used. Note that the PSD Sx(f) is centered about the carrier frequency with a 
spectral width of 2 f m = ~ . 

Assuming the scattering distribution of equation 5.16, it can be shown that the n th spectral moment of 
x(t) is given by [1] 



a n = a 0 (2n f m ) n q n {x) (5.18) 

where a 0 = \a 2 and 

q„(x)= 77 . [ e xcos0 cos " 9 d9 (5.19) 

nlo(x) Jo 

Specifically, based on Equation 5.18 and Equation 5.19, we obtain 

ai = a 0 {2jr (5.20a) 

kix) 

(0 ( \2 ( !o(x) + h(x)\ 

a 2 = a 0 {2n f n ) — . (5.20b) 

V 2I 0 (x) ) 

In the case of isotropic scattering (x = 0), the first two spectral moments reduce to 

= 0 and a 2 = 2 a 0 (jtf m ) 2 (5.21) 

These statistics will be used further in the design of velocity estimators and in the evaluation of their 
performance. 
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5.3 Principles of Mobile Velocity Estimation 

We can see from Equation 5.1 and Equation 5.9 that the information on mobile velocity is contained in 
both the envelope and phase of the received signal. The abundance of existing velocity estimators is a 
reflection of the numerous ways of extracting this information. The velocity of a mobile unit is generally 
estimated by exploiting the statistics of the envelope, phase, in-phase, and quadrature phase components 
or any other feature of the received signal. The approach for deriving a velocity estimator is based on 
the principle that there exists an exact, nonrandom relationship between the mobile velocity and the 
statistics of the Rician component s (t) in the presence of isotropic scattering. Several exact and equivalent 
expressions of the velocity can be derived, depending on which feature of s(f) is being considered, i.e., 
envelope, phase, in-phase, or quadrature phase, and which order of the statistic is chosen, i.e., first order, 
second order, etc. Below are some examples of these exact expressions of velocity. 



5.3.1 Examples of Derivations of the Velocity 

Below are some examples of how the mobile velocity can be calculated from the statistics of various features 
ofs(f). 



Example 5.1 Level Crossing Rate Method 

The envelope phase description of the Rician fading s(t) is 

s(t) = r(t) cos (2nf c t + i/r(t)) (5.22) 



where r(f) and i/r (t) represent the envelope and phase ofs(f), respectively. 

The level crossing rate (LCR) is defined as the average number of upcrossings per second the envelope 
r(t) makes of a predetermined level R„. It is obtained as [43, p. 77] 



,oo I / \ 

LCR r (R 0 ) = rp(r = R 0 ,r)dr = — W exp (— -3- ) (5.23) 

Jo do \ 2n \ 2a 0 J 

where p(r,r) denotes the joint probability density function (pdf) of r(t) and its time derivative f(f), 
a 0 = f er 2 , and a 2 is the second spectral moment of x(t), given in Equation 5.20b. If the level R 0 is chosen 
equal to s/2a 0 , the LCR of level R 0 becomes 



LCR r (Ro) 




(5.24) 



In the case of isotropic scattering, a 2 = 2ao(nf m ) 2 (Equation 5.21). It follows that Equation 5.24 reduces to 



LCR r (Ro) = \J~2jt- e 1 
X 



(5.25) 



The mobile velocity can then be exactly expressed in terms of the average number of upcrossings the 
envelope r(t) makes of level R 0 = s/2ao as 

ke 

v=—=LCR r (R 0 ) (5.26) 

spin 



Example 5.2 Rate of Maxima Method 

The rate of maxima (ROM) of a given process r (f ) is defined as the average number of maxima per second 
of r(f). It can be derived as [53] 



ROM r = 




0 ,r)df 



(5.27) 



where p(f,r) denotes the joint (pdf) of the first and second derivatives of the process r(t). 
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Using the results in [53] , the ROM of the envelope r ( t ) of signal s (f ) in Equation 5.22 can be obtained as 



ROM r = 




(5.28) 



which, in the case of isotropic scattering, and following the approach of Example 1, reduces to 



ROM r = 1.5651 — 
X 



(5.29) 



From Equation 5.29, the mobile velocity can then be exactly expressed in terms of the rate of maxima of 
the envelope of s (f) as 



v = 0.6389 XROM r 



(5.30) 



Example 5.3 Zero Crossing Rate Method 

The zero crossing rate (ZCR) of a given process a(t ) is defined as the average number of positive-going 
zero crossings per second of that process. It is given by 

poo 

ZCR a = / ap(a = 0,a)da (5.31) 

Jo 

In [54] , the ZCR definition was applied to the in-phase component of signal x(t) defined in Equation 5.1 1 
and Equation 5.15 and was shown to be equal to 



ZCR Xi = 




which in the case of isotropic scattering reduces to 



ZCR Xi = 



v 

sfix 



The mobile velocity can then be obtained using the following simple expression: 

v = \fl\ZCR Xi 



(5.32) 



(5.33) 



(5.34) 



The same result can be obtained using the quadrature component of x(t) instead of the in-phase com- 
ponent. 



Example 5.4 Covariance-Based Method 

For the covariance (COV)-based method, consider a process ri(f) = r 2 (f) where r(t) is the envelope of 
the Rician fading s(t). We can show that for a given time lag r 

E{(n(t+r) - n(t)) 2 } = 2R n (0) - 2R ri (r) (5.35) 



where R r ,(r) = £{ri(f + r)ri(f)} is the autocorrelation function ofprocess ri(f). 
In [7], it was proved that in the presence of isotropic scattering 



£{(u(f + r) - ri(f)) 2 } 



RrA 0) 




2 



(5.36) 



where x t is the sample spacing in seconds per sample and R ri (0) is the variance of the squared envelope. 
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The mobile velocity can therefore be extracted from Equation 5.36 as 



= 1 2R ri (0) - 2RJ^) 

2nr t \j R n { 0) 

where R r , (r) is the autocorrelation of the envelope square |s (f)| 2 for a given time lag r. 

Example 5.5 Instantaneous Frequency-based Method 

The IF of s (f ) is defined as [14] 



(5.37) 



ft At) 



1 rfi/r(f) 
2jt dt 



where ^r(f) is the phase of s(f). In Appendix A, we show that the first moment of |/i iS ] is 



(5.38) 



v / K\ -k 

(539) 

where | | stands for the absolute value operator, /) jS is the IF of the multipath component s (f ), and K is the 
Rice factor. The velocity of a mobile unit can therefore be expressed in terms of the first-order moment of 
the absolute value of the IF of the Rician fading component s (f) as 

v = V2XI~ 1 e*E{\fiA} (5.40) 



5.3.2 Examples of Velocity Estimators 

The expressions of the velocity in all the above examples are exactly equivalent. If we were to calculate 
the velocity using the LCR of the envelope of s(f) (Equation 5.26) or the first-order moment of the IF of 
s ( t ) (Equation 5.40), we would obtain exactly the same result. These exact expressions become estimators 
for two main reasons. The most important reason is that, in reality, we do not have access to s(t), and 
the scattering distribution is not necessarily isotropic. The exact expressions obtained from exploiting the 
statistics of features of s(t) in an isotropic environment become approximate in a nonisotropic environ- 
ment, and in the presence of shadowing and additive noise. The second reason is that some parameters 
required to compute the exact expression of v must also be approximated. For example, in Equation 5.40, 
since exact values for E { | /) jS | } or AT are not available, estimated values must be used, resulting in another 
source of error. 

In the presence of shadowing, additive noise, and nonisotropic scattering, the previous exact expressions 
become the following estimators. 

Example 5.1 LCR Estimator 

From the expression of the velocity in Equation 5.26, we can derive an estimator 

v LC r= ~^=LCR r (R 0 ) (5.41) 

s/2jt 

where LCR r (Ro) represents the number of level of upcrossings of the envelope of y(t) and where Rq = 
x/ 2 fl 0 ) with 

do = y J yf(t)dt — J yi(t)dt'j (5.42) 

where is the in-phase component of the received signal. 
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Example 5.2 ROM Estimator 

From the expression of the velocity in Equation 5.30, we can derive an estimator 

Vrom = 0.6389 XROM r (5.43) 

where ROM r is the rate of maxima of the envelope of the received signal y(f). 



Example 5.3 ZCR Estimator 

From the expression of the velocity in Equation 5.34, we can derive the following estimator: 

Vzcr = V2XZCR Zi (5.44) 

where ZCR Zj is the zero crossing rate of z;(f) = y;(f) — y, (f), where y, (f) is the in-phase component of 
the received signal y(f) and y; (f) an estimate of its mean. 



Example 5.4 COV-Based Estimator 

From the expression of the velocity in Equation 5.37, we can derive the following estimator: 



Vcov 



X I V 
2ttt, y R( 0) 



(5.45) 



where r f is the sample spacing in seconds per sample, R( 0) is the variance of the envelope square of the 
received signal, and 



V = 2(R(0) — R(t)) 



(5.46) 



Example 5.5 IF-Based Estimator 

The IF-based estimator is obtained from Equation 5.40 as 

v IF = XV 2 I 0 1 (jj ei < \f Uy (t)\ > (5.47) 

where 

< \fi,y(*)\ >= j J I fi, y (t)\dt 

is the time average of the IF estimate of the received signal over the duration T of the recorded signal, and 
K is the estimated Ricean factor. 

All the above-mentioned and other existing estimators exhibit different performances in the presence 
of additive noise, shadowing, and nonisotropic scattering. For example, in the presence of shadowing, 
the envelope of s (t) is distorted while its IF is not. As a consequence, any estimator based on the IF of the 
received signal would be expected to be more robust to shadowing than estimators based on the envelope 
of the received signal. The approach for evaluating the performance of velocity estimators is decribed in 
the next section. 



5.4 Performance Analysis of Velocity Estimators 

The examples of the previous section showed that several exact expressions of the velocity of a mobile unit 
can be obtained from exploiting different statistics and different features of the Rician component s(f). 
In the presence of additive noise, shadowing, and a nonisotropic environment, these exact expressions 
are used as velocity estimators with different properties and performances. Any performance analysis 
of a velocity estimator should evaluate its behavior in the presence of shadowing, additive noise, and 
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nonisotropic scattering. This approach is illustrated below for the IF-based estimator but can be applied 
to any other estimator. 

5.4.1 Effect of Shadowing 

In this section, the effect of shadowing is studied in the absence of additive noise. The model for the 
received signal is therefore y(t) = m(t)s (f). We observe that m(t) causes a distortion in the amplitude and 
quadrature components of the received signal. Thus, the performance of all velocity estimators that are 
based on the statistics of the envelope and quadrature components of the received signal deteriorates in 
the presence of shadowing. However, since the presence of shadowing changes only the amplitude of the 
received signal [43, p. 203; 62, p. 91], it produces no phase distortion. Hence, the IF of the received signal is 
not affected by m(f);i.e., both signals s(f) = r(f) cos(2jr/ c f+i/r(f))andy(f) = m(t)r (t) cos(2 jt f c t+i/s(t)) 
have the same IF, 4- x/r(t ). The IF-based estimator is therefore robust to shadowing and to any other 
amplitude distortion, such as path loss. 



5.4.2 Effect of AWGN and Nonisotropic Scattering 

5.4.2. 1 Derivation of the Normalized Bias 

In this section, the performance of the IF-based estimator is investigated in the presence of additive 
white Gaussian noise (AWGN), nonisotropic scattering, and the absence of shadowing. The performance 
criterion is the normalized bias of the estimator. Ideally, the variance of the estimator should also be 
derived. However, at the time of publication, there is no closed-form expression of the variance. This is 
also the case for all existing estimators. The variance will be investigated in the next section using Monte 
Carlo simulations. 

The normalized bias of estimator v IP in Equation 5.47 is defined as 

Vrc 

e(x,y,K)= — ~ 1 
v 

J2X , / K\ k , 

= —I 0 - 1 l T )eiE{\f i , y \}-l (5.48) 

where y is the signal-to-noise ratio (SNR) and 

E{\f, y \} = t- / W P^)dif (5.49) 

2lt J- oo 

p(iir) is given by [54] 
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Q 2 

P = tt 
2b 0 

B = b 0 b 2 - b\ 

b 2 - 2biir + b(,\fr 2 

U ~ 2 B 

_ Q 2 (b 2 - b i^) 2 
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(5.51) 



(5.52a) 

(5.52b) 

(5.52c) 
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and b„ is the n th spectral moment of ya(t) = m(t)x(t) + n(t). yo(t) is the received signal in the absence 
of an LOS component. 

In the absence of shadowing, m{t) is a constant M. Since both x{t) andn(f) are zero-mean, independent 
processes, the n th spectral moment of yo(t) can be written as 

POO 

b n = ( 2nT / (/- f c nM 2 S x (f) + S N (f))df 
Jo 

= M 2 fl„ + (l-(-l)" +1 ) NoJT B' 1+1 , (5.53) 

4(n + 1) 

n > 0 



5.4. 2. 2 Performance in the Presence of AWGN and Isotropic Scattering 

In the case of isotropic scattering (/ = 0, Figure 5.5(a)), E {| f uy \ } is given by (Equation 5.A5) 






2 7T V ^0 

The normalized bias of Equation 5.48 can then be written as 

X 



e(0 ,y,K) = 






K \ K- p 

y' eT -i 



(5.54) 



(5.55) 



sfinv 

where p is defined in Equation 5.51. After using Equation 5.18 and Equation 5.53, Equation 5.55 becomes 



e(0, y, K) 



y + i 



Ky 

2(y + 1) 




K 

e 2(y+l) — 1 



(5.56) 



Figure 5.6 and Figure 5.7 show the effect of AWGN on the IF-based estimator with respect to y and v, 
for two different values of K. The bandwidth B 0 is assumed to be B 0 = = 170 Hz, which allows 




FIGURE 5.6 The normalized bias of the IF-based velocity estimator as a function of SNR for three different MS 
velocities, assuming K = 0. 
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FIGURE 5.7 The normalized bias of the IF-based velocity estimator as a function of SNR for three different MS 
velocities, assuming K = 10. 

for velocities up to v max = 100 km/h at f c = 900 MHz. We observe that for low velocities and low SNR, 
the bias is significant. The normalized bias is also seen to increase slightly as K increases. This is also the 
case for some existing estimators [7], However, Appendix C shows that the IF estimator allows for an 
improvement in SNR by a factor of 3 In where /i is the lower cutoff frequency of the output low-pass 
filter at the output of the IF estimator. For typical values, B 0 = 170 Hz and /) = 10 Hz, the improvement 
is about 8 dB. This means that the normalized bias of the IF-based estimator for an input SNR of 10 dB 
is the same as that of the ZCR method for an input SNR of 18 dB. Thus, even for small velocities, the 
IF-based estimator will still exhibit a negligible bias at low SNRs due to the improvement of SNR at the 
output of the IF estimator. 

5.4. 2.3 Performance in the Presence of Nonisotropic Scattering 
and in the Absence of AWGN 

The effect of the scattering distribution is investigated here in the absence of noise (y — > oo). In this case, 
using Equation 5.50, the normalized bias in Equation 5.48 in the absence of additive noise, e(x, oo, K), 
can be computed numerically. Figure 5.8 shows the effect of the scattering distribution on the IF-based 
estimator for different values of / . Because the estimator in Equation 5.47 was derived in the case of isotropic 
scattering, the error in estimating the MS velocity increases as / increases. Comparing Figure 5.8(a) and 
(b), we observe that there is a negligible increase in normalized bias when K increases. This suggests that 
e(x,oo, K) can be approximated by e(x,oo,0). Using Equation 5.48 and Equation 5.A5, we obtain 

e(/,oo,K) ~ e(x,oo,0) 

= \/2 qiix) — 1 ( 5 . 57 ) 

Appendix B shows that the ZCR-based velocity estimator (Equation 5.44) has the same normalized bias 
as the IF-based estimator (Equation 5.47) in the presence of nonisotropic scattering. The ZCR-based 
velocity estimator is generally more robust than the LCR and covariance-based methods in the presence 
of nonisotropic scattering [7]. It then follows that the IF-based velocity estimator also outperforms the 
above-mentioned estimators in the presence of nonisotropic scattering. 
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FIGURE5.8 ThenormalizedbiasofthelF-basedvelocityestimatorasafunctionof/: (a) forK = Oandv = 40 km/h, 
(b) for K = 10 and v = 40 km/h. 

In summary, unlike the ZCR, LCR, and COV-based methods, or any method based on the envelope, in- 
phase, or quadrature phase of the received signal, the IF-based velocity estimator is robust to shadowing. It 
also exhibits superior performance in the presence of AWGN due to the improvement in SNR introduced 
by the IF estimator. In the presence of nonisotropic scattering, the IF-based estimator has the same 
performance as the ZCR and is generally more robust than LCR and COV-based methods. The effects 
of additive noise, nonisotropic scattering, and shadowing combine in a nonlinear way and are easier to 
evaluate using simulations. This is the object of the next section. 
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5.5 Performance Analysis Using Simulations 



In Section 5.4, the performance of the IF-based estimator was evaluated analytically in the case of re- 
strictive assumptions. For example, the IF-based estimator was evaluated in terms of the bias introduced 
in the presence of additive noise, shadowing, and nonisotropic scattering considered separately. It is also 
important to investigate the performance of an estimator when all these effects are present at once, which 
is the case in real life. Analytical derivations can, in that case, become intractable, and it is often easier to 
use simulations. Also, the IF-based estimator was evaluated in terms of bias only. However, the variance 
of an estimator is another important measure of performance. At time of publication, in the case of the 
IF-based estimator or any other velocity estimator, there is no closed-form expression for the variance. 
In this case, Monte Carlo simulations are used. Finally, the analytical derivations of Section 5.4 assumed 
perfect knowledge of the first-order moment of the IF of the received signal and of the Rice factor. A more 
accurate performance study would require taking into account the statistics of the estimators used for the 
first order of the IF and for the Rice factor. However, the derivations would again become too complex 
and simulations are a better alternative. 



5.5.1 Simulations of the Received Signal 

Using Jakes method for simulating the multipath component of the received signal, the in-phase and 
quadrature components of x(t) defined in Equation 5.15 can be modeled as [34] 



x, ( t ) = ^2 cos ( 27 r f m t cos f ) + 4y 



n = 1 
N 



;(f) = 5Z Sill ( 2jr /™ tC °s(“^"“ ) + 



(5.58a) 

(5.58b) 



where N is the number of waves arriving at the MS antenna and (« = 1, 2, . . . , N) are assumed to be 
i.i.d. uniformly on (0, 2 jt). The process L ( t ) of Equation 5.2 is modeled as a zero-mean Gaussian process 
with PSD Sl(v) given by 



S L (v) 



2d 0 crl 

1 + (2nvd 0 ) 2 



(5.59) 



where d 0 and ay are the correlation length and variance of L , respectively. Let v max denote the maximum 
spatial frequency of Si (v), and D = vT, with T the total duration of the simulated signal. A model for L 
is [49] 



/-i r 
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(5.60) 



where 9i(k = — /, — / + 1, ...,/) are assumed to be i.i.d. uniformly on (0, 2j r), / = Dv max , and 



C = 



1 

Dal 



7-1 




k + 1/2 \ 
D 



(5.61) 



In these simulations, the carrier wavelength is A = 1/3 (f c =900 MHz), and the exponent of the dis- 
tance dependence is a = 2. The standard deviation of lognormal shadowing ranges from 6.5 to 8.2 dB 
at 900 MHz in urban areas for microcells [47]. Here, the standard deviation and correlation length of 
the lognormal shadowing are ay = 8 dB and rf 0 = 50 m, respectively. Parameter P 0 is assumed equal 
to 1 W. 
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FIGURE 5.9 Block diagram of the IF-based MS velocity estimator. 



5.5.2 Simulation Results 

We consider that only a finite duration of the received signal is observed, and that the IF and moments 
of the signal are estimated using a finite duration of the received signal. The performance of the IF-based 
estimator is then compared with that of other estimators. 

For each value of the MS velocity, the estimator P IF in Equation 5.47, is computed by using the approach 
illustrated in Figure 5.9. A limiter is used to remove spurious amplitude variations from the received 
signal without destroying the information in the phase. It is followed by a bandpass filter of suitable 
bandwidth B 0 , in order to extract the constant-amplitude signal. A value of B 0 = 170 FIz is chosen for 
the system bandwidth. This process is widely used in FM receivers [18, pp. 205-207], The signal z(f) is a 
constant-amplitude frequency-modulated signal and can be expressed as 

z(t) = Acos(2nfot + <p(t)) + n 0 (t) (5.62) 



where A is the maximum amplitude at the limiter output and n 0 ( t ) is the bandpass noise at the output of 
the filter. The IF of the signal z(f) can be estimated using [15] 



1 zMz' q {t) - z](f)z,(f) 
kz(t ~ 2 jt z?(f) + zg(f) 



(5.63) 



with z;(f) = Acos(</>(f)) + n„i(t) and z q (t) = Asin(0(f)) + n oq (t). In practice, however, this method is 
difficult to implement as coherent detectors are needed to recover the in-phase and quadrature components 
of z(f). In order to overcome this problem, a balanced discriminator can be deployed to estimate the IF of 
z(f) [18, p. 216]. This is the method used in these simulations. The samples of the signal z(f) are passed 
through an finite impulse response (FIR) differentiator. The envelope of the output of the differentiator is 
proportional to Therefore, unlike other estimators such as the ZCR estimator, the IF-based estimator 

does not require coherent detectors. The Ricean K factor is estimated using the estimator of Equation 
5.70 derived in the next section. The IF-based estimator of Equation 5.47 can then be calculated using the 
envelope of the estimated IF and the estimated K factor. 

The performance of the IF-based estimator in the presence of shadowing and additive noise was com- 
pared with that of the ZCR, LCR, ROM, and COV-based methods by computing the average normalized 
bias and variance over 100 realizations of the received signal. The sensitivity of the COV-based method 
to additive noise can be reduced significantly by choosing a large time lag. However, a large time lag 
reduces the accuracy of the estimator itself. Here, the time lag in the COV-based method was chosen to 
be 2.5 msec, which according to our simulations gives an accurate estimate of the MS velocity in the ideal 
case. In order to limit the delay in obtaining velocity estimates for real-time implementation, we used a 
window of 1 sec of the simulated signal to estimate the unknown velocity. As the window length increases, 
the bias introduced by estimating E { | f ) iy \ } using the time average of the estimated IF decreases, and the 
delay in obtaining velocity estimates increases. 
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FIGURE 5.10 Normalized bias as a function of the MS velocity in the presence of shadowing. 

Figure 5.10 shows the effect of shadowing on each of the velocity estimators in the absence of additive 
noise and in the case of isotropic scattering. We observe that in the presence of shadowing, the normalized 
bias in estimated velocity for the IF-based estimator is negligible compared to the other estimators. This 
confirms that the IF-based estimator is robust to shadowing. 

Figure 5.11 shows the effect of additive noise in the absence of shadowing, assuming isotropic scattering. 
We observe that the IF-based estimator is generally more robust than the other estimators. It exhibits a 




FIGURE 5.1 1 Normalized bias as a function of the MS velocity in the presence of AWGN for SNR =10 dB. 
Copyright © 2005 by CRC Press LLC 







FIGURE 5.12 Normalized bias as a function of the MS velocity in the presence of AWGN, shadowing, and isotroic 
scattering for SNR = 10 dB. 



smaller bias than the other estimators for small velocities. This is due to the improvement in SNR at the 
output of the IF estimator (see Appendix C) . As the MS velocity increases, the system bandwidth approaches 
the actual Doppler shift, and therefore, the difference between the normalized biases of all the estimators 
becomes particularly small. Finally, Figure 5.12 to Figure 5.14 represent the respective normalized bias, 




FIGURE 5.13 Variance as a function of the MS velocity in the presence of AWGN, shadowing, and isotroic scattering 
for SNR = 10 dB. 
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FIGURE 5.14 Mean square error as a function of the MS velocity in the presence of AWGN, shadowing, and isotropic 
scattering for SNR =10 dB. 



the variance, and the mean square error of the above-considered estimators in the presence of additive 
noise, shadowing, and isotropic scattering. They indicate that the ROM estimator exhibits slightly lower 
variance than the IF-based estimator, but higher bias and mean square error. 



5.6 Rice Factor Estimation 



5.6.1 Existing Methods 

The value of K influences most local power and mobile velocity estimators [7]. It is also a measure of the 
severity of the fading, K = 0 being the most severe Rayleigh fading and K =oo indicating that there is no 
fading. Hence, knowledge of the Ricean K factor is a good indicator of the channel quality [25]. 

Traditional methods for estimating the Ricean K factor are based on measurements of the received 
power. In [25], the probability distribution of the measured data is computed and compared to a set of 
hypothesis distributions using a suitable goodness-of-fit test. In [46], a maximum likelihood estimate of 
K is obtained using an expectation-maximization (EM) algorithm. Both methods have a high degree 
of complexity, which makes them relatively time-consuming [24] . A simple and rapid method for K 
estimation uses two estimated moments of the received power [24]. The same estimator has been derived 
in [66] using an approach based on the covariance of the received power. In [63], the Ricean K factor is 
estimated using the first two moments of the envelope of the received signal. However, this method did not 
receive much attention, as it did not provide a closed-form expression for K . However, Tepedelenlioglu et al. 
[64] recently proposed a general class of moment-based estimators for K, which included the estimator of 
[63] . They derived the asymptotic variance of each member of the family and showed that the asymptotic 
variance of the estimator in [63] is the closest to the Cramer-Rao lower bound. In this section, we present 
two explicit and simple estimators for K based on an approximation of the method in [63] . We also show, 
using simulations, that one of the proposed estimators has a lower mean square error than the estimator 
in [24] and [66]. 
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5.6.2 Envelope-Based Estimators 

The ratio of the first two moments of the envelope, r(f), of the Rician fading is [53] 



a E{r} 



2ym 



(n-ni„(f) +*/,(§) 



(5.64) 



Using Equation 5.64, K can be obtained by first estimating the first two moments of r(f), and then 
solving Equation 5.64 numerically for K . This was the method originally proposed in [63 ] . Elowever, since 
Equation 5.64 does not provide an explicit formula for K, it has not been used as an estimator. 

To overcome the aforementioned drawback, we approximate the right side of Equation 5.64, and derive 
two explicit estimators for K. We first rewrite Equation 5.64 as 



g(K) 
r K + 1 



(5.65) 



where g ( K ) is defined as 



g (K)=^e~^VK + 1 



u + io/oifUWf 



(5.66) 



Figure 5.15 plots g(K) for K e [0, 100]. It suggests that g(K) can be well approximated with a low-order 
polynomial function gN(K), where 



N 

g N (K) = Y Pm k ' {5.67) 

>=o 

The coefficients (z = 0, 1, . . . , N) are computed by fitting gn(K) to g(K) in a least squares sense. It 
can be shown that for N > 2, p; N , i > 2 is of order 10 -5 . We choose to approximate g(K) with a linear 
and quadratic function, gi(K) and giiK), respectively. 




FIGURE 5.15 g(K) in Equation 5.66 vs. K. We observe that g{K) can be well approximated by a low-order 
polynomial. 
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(5.68) 



Replacing g(k) by its linear approximation in Equation 5.65 leads to the following estimator: 

pp Er - poi 

Jti = ~ 

Pn - E r 

In Equation 5.68, E r is the ratio of the estimated first and second moments of the envelope of the received 
signal, i.e., E r = R/\fs 1 where, assuming local ergodicity, 

R" = ^ J R"(t)dt, n = 1,2 

and T is the duration of the observed signal. By fitting g i ( K ) to g ( K ) in a least squares sense, we find 

poi = 0.7967, p n = 0.9969 (5.69) 

Following the same approach, g(k) is replaced by its quadratic approximation in Equation 5.65, which 
leads to the following estimator: 



K 7 = 



E r — pn + V ( E r — pn) 2 + 4 p 2 i(E r — poi) 
2p22 



(5.70) 



where 



P02 — 0.8293, pn — 0.9866, p 22 — 0.0005 
Parameter E r in Equation 5.64 can be approximated by 

Poi + PnK 



E r ~ Er I = 



or 



E r ~ Er? = 



K + 1 

P 02 + PnK + P 22 K 2 
K + 1 



(5.71) 

(5.72) 

(5.73) 



depending on the degree of accuracy required. 

Figure 5.16 plots Equation 5.64 and its approximations (Equation 5.72 and Equation 5.73). It verifies 
that E r 2 is a better approximation than E r \. 




FIGURE 5.16 E r vs. K (Equation 5.64) and its approximations (Equation 5.72 and Equation 5.73). 
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FIGURE 5.17 ( Continued ). 



5.6.3 Simulation Results 

In this section, the performance of the proposed estimators is compared to that of the COV-based estimator 
[24] , [66] . The performance criterion is the mean square error of the estimators. The COV-based estimator 
is defined as [66] 



Kcov 



a y ~ %l 2 (°) + °W a y ~ %l 2 ( 0) 
C |y|2 (0) 



(5.74) 



where C| y | 2 (r) refers to the covariance of|y(f)| 2 and er 2 is the power of the received signal. 

The mean square error (MSE) ofthe three estimators was computed over 1000 realizations of the received 
signal. Figure 5.17 shows the MSE of the estimators as a function of SNR. We observe that for SNRs larger 
than 15 dB, the estimator given by Equation 5.70 is superior to the estimators given by Equation 5.68 and 
Equation 5.74. 



5.7 Application on Handover Performance 



5.7.1 Handover Decision Algorithms 

In cellular, microcellular, and picocellular systems, as a vehicule crosses the cell boundary between two BSs, 
control has to be transferred from the current BS to the target BS. This process is referred to as handover 
or handoff and is generally categorized into soft handoff and hard handoff processes [70, 75]. In hard 
handover processes, only one BS is in control of an MS at a given time. Soft handover processes require two 
or more BSs to be allocated to a single MS, in order to prevent the received signal power from dropping 
below a given threshold level. Soft handover processes are currently specific to CDMA systems where 
power control is critical [75]. The aim of this section is to study the effect of a biased velocity estimation 
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on the performance of handover decision algorithms. The derivations are done in the case of hard handoff 
processes but can also be adapted to soft handoff processes. 

Consider an MS traveling from its current base station, BS1, to a target base station, BS2 (Figure 5.18). 
The decision to hand off is based on the received signal strength from both base stations. In the absence 
of shadowing, fading, and neglecting the effect of additive noise, the received signal envelopes are only 
affected by path loss. The path loss is inversely proportional to d a , where d is the distance to the BS and a is 
a parameter determined by the environment. At mid-distance from the two BSs, the received strengths are 
identical, or equivalently, the difference between the two received strengths is zero. A decision to hand off 
from the current BS to the target BS could then be made when the sign of the difference between the two 
received strengths changes. However, in the presence of shadowing or fading, the sign can change several 
times around the cell boundary, which could result in a large number of unnecessary handovers. This is 
called the ping-pong effect. The effects of fast fading are significantly reduced through time-averaging the 
received signals over a window of size L. In [62], it was recommended that the window size corresponds 
to a distance of 40A, hence to a period of time T = 40 k/v. Assuming the signal strengths are measured 
every T s seconds, the time averaging should occur over L = [^] samples, where [x] < x < [ x ] + 1. 
The window size is therefore velocity dependent and is smaller for higher velocities than for smaller 
velocities. 



5.7.2 Effect of an Error in Velocity Estimation on the System's 
Quality of Service 

The quality of service of a system requires that the probability of a call being dropped must be less than 
pout) where p ollt is referred to as outage probability. Here we investigate how an error in velocity estimation 
Sv would impact the probability of lost calls. 

Every sampling time T s , the MS of Figure 5.18 records measurements of the received signal strengths 
from BS1 and BS2 and computes a time average of the two received signal envelopes based on current and 
past measurements. Let A y\ (k) and A y 2 (k) be the k th sample of the time-averaged signals from BS1 and 
BS2, respectively. Assuming a measurement time k > L, Ayi(k) and A y 2 (k) can be expressed as 

, k 

A yi(k)=— ^2 y> 0’)> * = { 1.2} (5.75) 

j=k—L + 1 

where y,(j) is the signal strength received from BSi at time j . 

An error, Sv, in velocity estimation leads to the choice of a temporal window 



L' = 



40A 

(v + Sv)T s 



= L+SL 



(5.76) 
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FIGURE 5.18 We consider an MS traveling from its current base station, BSI, toward a target base station, BS2. 
Copyright © 2005 by CRC Press LLC 




where SL is a positive integer when Sv is negative and a negative integer when <5v is positive. It follows that 
the averaged signal strengths in Equation 5.75 become 



A >i(fc) = 



1 

L+SL 



y yi(j), i = {1,2} 



j=k-L-8L + 1 

Assuming an error 8L such that 8L/L « 1, Equation 5.77 can be written as 
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(5.77) 
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where Ay,- ( k ), i = {1,2}, is defined in Equation 5.75. An error SL in window size leads to an error <5e,- in 
average signal strength, where <5e,- is defined as 
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where i = {1,2}. 

It is assumed, in handover decision algorithm design, that the process of time averaging has smoothed 
the effect of fast fading but not that of shadowing [ 5 1 ] . As a consequence, merely handing over from BS 1 
toward BS2 and back to BS1 when the difference Ayi(fc) — A y 2 (k) changes sign would still result in a high 
number of unnecessary handovers. A more appropriate handover decision algorithm consists of handing 
over from BS1 toward BS2 when A yi(k) — A y 2 {k) < —h, and back to BS1 when A yi(k) — A y 2 (k) > h, 
where h is a hysteresis value determined by the system and by the environment (see Figure 5.19) [70]. The 
choice of h can be critical in terms of minimizing two conflicting criteria. If h is too small, shadowing is 
still dominant, resulting in a high number of unnecessary handovers. If h is too large, the handover process 
is delayed and a call could be lost. 

A call is lost when the base station is BS1 and A yi(k) drops below a given threshold th determined 
by the sensitivity of the receiver, or similarly, when the base station is BS2 and A y 2 (k) is below th. The 
probability of a call being dropped at time k is then [70] 

pLc(k) = prob(BSl(k)fkAyi(k) < th) + prob(BS2(k)tkAy 2 (k) < th) (5.87) 

If an error Sv occurs in the estimation of the velocity, a decision to hand off from BS1 toward BS2 
would now be made when A 'yi{k) — A'y 2 (h) < —h, and back to BS1 when A r yi(k) — A 'y 2 (k) > h. 
Equivalently, a hand off toward BS2 would occur when A yi(k) — A y 2 (k) < —h + Sh, and back to BS1 
when Ayi(fc) — A y 2 {k) > h + Sh, where Sh = 5ei(fc) — Se 2 (k). This could be interpreted as the decision 
to hand off in a system where the mobile velocity is known or is estimated in a very accurate way, but 
Copyright © 2005 by CRC Press LLC 



BS1 to BS2 



BS1 to BS2 



FIGURE5.19 A handover is processed toward BS2 when Ayi(k) — Ay 2 (k) < —h, and back to BS1 when Ay! (k) — 
Ay 2 (k) > h, where ?! is a hysteresis value determined by the system and by the environment. 

where two distinct hysteresis values, h\ = — h + Sh and h 2 = h + Sh, are used for handing over from 
BS1 to BS2 and from BS2 to BS1, respectively, where h i /it- It is assumed that the hysteresis value h 
was optimized so that the probability of lost calls resulting from a delay in handover is below the outage 
probability of the communications system. However, in the presence of an error in velocity estimation, if 
Sh is negative, a delay in handing over from BS1 toBS2 occurs, and USh is positive, a delay in handing over 
back to BS1 occurs. Either way, a delay in handover occurs, and depending on the value of Sh, this delay 
could contribute to increasing the probability of losing a call. There is another way an error in velocity 
estimation can contribute to an increase in the probability of lost calls. Because of the error in velocity 
estimation, the probability of losing a call at time k is now 

pLc(k ) = profr(BSl(/0&AyiW < th — Sei(k)) + prob{BS2{k)&.Ay 2 {k) < th — Se 2 {k)) (5.88) 

Equation 5.88 can be interpreted as the probability of lost calls in a system where the mobile velocity is 
known or is estimated in a very accurate way, but where the sensitivity of the receiver is different. If S ( k ) 
and Se 2 (k) are positive, the receiver can be considered as having better sensitivity. However, if either Sei(k) 
or Se 2 (k) is negative, the sensitivity of the receiver could be considered worse, which could contribute to an 
increase in the probability of lost calls. In [29], it was found that an increase of 20% in velocity estimation 
error results in the probability of lost calls increasing by 50%. 

5.8 Conclusions and Perspectives 

This chapter outlined the principles and methods of velocity estimation, provided examples of velocity 
estimators, and gave a framework for analyzing the performance of velocity estimators. It also discussed 
the implications of erroneous velocity estimation on the wireless telecommunications systems and, more 
specifically, in the context of handover algorithm design. It was found that an error in velocity estimation 
could significantly deteriorate the quality of service of the system in terms of probability of dropped 
calls. The CDMA 2000 system is currently being implemented in Australia. Although the service provider, 
TELSTRA, sets the standards to be met by mobile phone manufacturers, it has no control over which or 
even if a velocity estimator is to be used by the telephone hand set manufacturer. As a result, two mobile 
phone sets from different manufacturers used at the same location can have different performances. This 
suggests the need, in the future, for more control over which standards are to be used by mobile phone 
manufacturers. 
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Appendix A: Derivation of Equation 5.39 

In this appendix, we derive the first moment of the envelope of | f us |. Using Equation 5.38, the expected 
value of | fj t y | can be written as 

1 I* OO j*7 X 7*00 poo 

E{\fi, y \}=— / / / / W P(r,\ls,f,ir) th]r dr d\[r dr (5.A1) 

^ Jo J —71 J — OO j —OO 



where p(r, \fr, r , xfr) is the joint pdf of r, i/r , f , x/r. In the case of isotropic scattering, the joint pdf is given 
by [54] 



pi(r,^,r, Jr) 



r 



2 



4n 2 b 0 b 2 




x exp ^ 


-1 

2b 0 


x exp ^ 


-1 

2b 2 



(r 2 — 2 Qr cos Jr + Q 2 ) 
(r 2 + r 2 Jr 2 )] 



(5.A2) 



Substituting Equation 5.A2 into Equation 5.A1, and since pi(r, Jr, r, Jr) is an even function of i Jr, 
Equation 5.A1 can now be written as 

j /*00 j*7t POO POO 

E{\fi, y \} = — / / / / Jr pi{r,Jr,r,Jr)dJr df dir dr (5.A3) 

^ Jo J —n J —oo Jo 



Using the following identities [23], 

/ OO 

exp 



-1 

2 V 2 



f 2 ] dr = \j2nb 2 



1 jr exp [ —j—r 2 Jr 2 ) djr = b 2 r 



exp I — Qr cos Jr) djr = 2i tI 0 ( J-r 



bo 



exp(-/r 2 x 2 ) I 0 (vx)dx = exp ( ) I 0 ( — -r 



2/x \8/r 2 



8 ff 2 



in Equation 5.A3, we obtain 



1 jb 2 
2jr V b 0 



£fl/i,rl} = — \/7f e xp( — 



(?)'«(?) 



(5.A4a) 

(5.A4b) 

(5.A4c) 

(5.A4d) 

(5.A5) 



where p = ^ . In the absence of noise ( y -> 00 ) and in the presence of isotropic scattering, bo and b 2 are 
equal to a 0 and 2jr 2 f 2 cio, respectively, and since K = ^4- , Equation 5.A5 becomes Equation 5.39. 



Appendix B: Effect of the Scattering Distribution 
on the ZCR Method 

In this appendix, we study the effect of nonisotropic scattering on the estimator in Equation 5.44, in the 
absence of shadowing and additive noise. In this case, ZCR ^ in Equation 5.44 is given by Equation 5.32. 
From Equation 5.53, in the absence of shadowing and additive noise, bo and b 2 become M 2 a Q and M 2 a 2 , 
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respectively. In this case, Equation 5.44 reduces to 



Nzcr = — \J qiix) (5.B1) 

It follows that the normalized bias for the ZCR-based velocity estimator is derived as 

e(x,oo, K) = \J 2q 2 (x ) ~ 1 (5.B2) 

which is the same as Equation 5.57. Equation 5.B2 proves that both ZCR and IF-based estimators have the 
same performances in the presence of nonisotropic scattering. 



Appendix C: SNR Improvement in the IF Estimator 

In this appendix, we prove that the IF estimator improves the SNR. If we assume that the power of the 
noise at the input of the IF estimator is N 0 B 0 , the SNR at the input of the IF estimator is derived as 



Vr = 




A 2 

2 N 0 Bq 



(5. Cl) 



Based on the results in [18, Section 10.3], when y R is more than a threshold level 3 , the signal y 0 „ f (f) at the 
output of the IF estimator, can be expressed as 



Youtit) = = /], s (f) +£(f) 



where § (f) represents the IF noise. The PSD of the IF noise is given by [ 18, p. 414] 



W) = 




where IK-) represents the rectangular pulse function. Thus, the power of the IF noise in the frequency 
band ( /i, ^ ) can then be computed as 



2 No ( Bq 3 



N ° = 2 1 1 s ^ {f)df = [t ~ ^ 



The power of /; (f ) in the frequency band ( , A ) can also be computed as [43, p. 256] 



1 \ _ r* B 2 B 0 



Therefore, assuming /) « B 0 , the SNR at the output of the IF estimator can be derived as 



Kn=(^) ^3y«ln^ 



(5.C2) 



It follows from Equation 5.C2 that the IF estimator improves the SNR by a factor 3 In 



3 This threshold is found to occur in the vicinity of 10-dB input SNR and varies in the approximate range of 6 to 13 
dB [67, p. 458]. 
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